function r = grad(h,k,x)
dx=4*h*x(1)*(x(1)^2+x(2)-11)+2*k*(x(1)+x(2)^2-7);
dy=2*h*(x(1)^2+x(2)-11)+4*k*x(2)*(x(1)+x(2)^2-7);
r = [dx, dy]';
end